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We analyze the properties of breathers (time periodic spa- 
tially localized solutions) on chains in the presence of alge- 
braically decaying interactions l/r s . We find that the spatial 
decay of a breather shows a crossover from exponential (short 
distances) to algebraic (large distances) decay. We calculate 
the crossover distance as a function of s and the energy of the 
breather. Next we show that the results on energy thresholds 
obtained for short range interactions remain valid for s > 3 
and that for s < 3 (anomalous dispersion at the band edge) 
nonzero thresholds occur for cases where the short range in- 
teraction system would yield zero threshold values. 
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The understanding of dynamical localization in classi- 
cal spatially extended and ordered systems experienced 
recent considerable progress. Specifically time-periodic 
and spatially localized solutions of the classical equations 
of motion exist, which are called (discrete) breathers, or 
intrinsic localized modes. The attribute discrete stands 
for the discreteness of the system, i.e. instead of field 
equations one typically considers the dynamics of de- 
grees of freedom ordered on a spatial lattice. The lattice 
Hamiltonian is invariant under discrete translations in 
space. The discreteness of the system produces a cutoff 
in the wavelength of extended states, and thus yields a 
finite upper bound on the spectrum of eigenfrequencies 
O g (phonon band) of small-amplitude plane waves (we 
assume that usually for small amplitudes the Hamilto- 
nian is in leading order a quadratic form of the degrees 
of freedom) . If now the equations of motion contain non- 
linear terms, the nonlinearity will in general allow to tune 
frequencies of periodic orbits outside of the phonon band, 
and if all multiples of a given frequency are outside the 
phonon band too, there seems to be no further barrier 
preventing spatial localization (for reviews see B] , || ) . 

To cope with breathers in lattice dynamics, one has 
to face the problem of (i) quantization of breathers, ii) 
breathers in the presence of acoustic phonon bands and 
iii) breathers in the presence of long range interactions 
(e.g. in ionic crystals). While (i) still lacks a full under- 
standing, the correspondence between classical breathers 
and quantum bound states is believed to be correct ( 
§, Hi, il , §, [§)• The case of acoustic breathers has 
been studied in one dimension in || and in two dimen- 
sions in JsJ , where it was shown that the resonance of a 
zero frequency component of the breather (static defor- 
mation) with the zero of the acoustic spectrum leads to 



an algebraically decaying lattice deformation, but not to 
a disappearance of the breather. 

As for the case of long range interactions, less results 
are known (but see e.g. jl(|, |[2]]). A general proof of 
existence of breathers in d-dimensional lattices with alge- 
braically decaying interactions was obtained in JTT[ with 
upper bounds for the spatial decay of the breather am- 
plitude. Namely for interactions decaying like l/r s with 
r being the distance from the breather center and s some 
power s > d, the breather amplitude is bounded from 
above by a power law a/r s with a being some nonzero 
constant^]. This result leaves us with two questions. 
First, how can one obtain contact with the case of short 
range interaction (basically s — > oo) where exponential 
localization takes place? And secondly, what is hap- 
pening to energy thresholds of discrete breathers in the 
presence of long range interactions? In the case of short 
range interaction, simple estimates of the far distance en- 
ergy of a breather solution yield the correct predictions 
for nonzero energy thresholds [jl4|. Obviously these es- 
timates would yield zero energy thresholds for all cases 
of long range interactions if the far distance energy is 
calculated with the help of a spatial decay a/r s for the 
breather (the far distance energy would be simply ~ a 2 , 
and would always tend to zero if a — ► 0, see also jf3|). 
We will resolve these puzzles in the following. 

At this stage it is appropriate to fix the class of Hamil- 
tonians to be considered further: 



H = J2 
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Here P; and Xi are canonically conjugated scalar mo- 
menta and displacements of a particle at lattice site I. 
The on site potential V(z) — Y^=2 ~jr z ^ can ^ e use d as a 
simple way to generate an optical phonon spectrum, and 
the interaction Wi{z) = J2^=2 ^jp-z 11 should incorporate 
longe range interactions with ^(l) — y~F- For small val- 
ues of Pi and Xi the classical Hamiltonian equations of 
motion Xi = , Pi = — can ^ e linearized in Xi. 
The corresponding eigenvalue problem when solved for 
plane waves Xi(t) ~ exp 1 ^' - ^'- 1 is given by 



1 This is in contrast to results obtained in jiJJ and JUsf , 
where for d — 1 and s > 3 exponential decay was obtained. 
Actually both decay laws hold, see below. 
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n 2 q = v 2 + 2C J2 — s C 1 - cos(qm)) . (2) 

m— 1 

Let us discuss the properties of E s {q) = > 0. First 
E s {q) is bounded from above for all s > 1 and periodic 
in q with period 2-7T. Most important is that E s (q) is a 
nonanalytic function in i.e. its n — (s — l)-st deriva- 
tive with respect to q is discontinuous at q = (when 
s is noninteger, (s — 1) < k < s). This follows already 
from the fact that the convergence radius of (^) is zero 
for nonzero imaginary components in q. Indeed for even 
integers s one finds (E s (q) — v 2 ) ~ B s (q/(2ir)) for 
< q < 27r. Here B s (z) is the Bernoulli polynomial 
of s-th order. Consequently at small q the expansion of 
E s (q) contains a term q s_1 which leads together with 
the periodicity of E s (q) to the mentioned nonanalyticity. 
For odd integers s the expansion of E s (q) contains a term 
q"^ 1 ln(q), and for noninteger s a term q s_1 (follows from 
d 2 E s (q)/dq 2 = -E s ^ 2 {q) + 2C((s - 2) with ((z) being 
the Riemann Zeta function) . Finally for small q the lead- 
ing term in the expansion of E s (q) is v 2 + C(^(s — 2)q 2 
for s > 3 and v 2 + 2Ca(s)q s ~ 1 for 1 < s < 3 with 
a(s) = (I — cos x) / x s dx [pL3[| . Note that the dispersion 
at the upper band edge (q — ir) is completely analyti- 
cal, and in leading order always proportional to (q — ir) 2 . 
Some of these results have been discussed at length in 
p6[ (see also original references therein). 

Now we can turn to the first problem of the spatial de- 
cay of a breather. In order to generate a breather solution 
we chose V4 ^ and all other anharmonic terms in V(z) 
and W{z) being zero. Since we can only simulate finite 
system sizes N , we use periodic boundary conditions. In 
that case we have to define a cutoff length in the inter- 
action which we chose to be N/2 (we will discuss the 
corresponding corrections to E s (q) later). We calculate 
breather solutions using a Newton algorithm (see (lTj , ||] 
for details). The results for s = 10,20,30 are shown in 
Fig.l. 
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FIG. 1. Breather solution at time t — with 

Pi(t = 0) = 0. The corresponding displacements (amplitudes) 
Xi(t — 0) are plotted versus lattice site. The nonzero model 
parameters are t>2 = V4 — 1, C = 0.01. The frequency of 
the solutions is chosen T — 4.7682. Circles: s = 10, squares: 
s = 20, diamonds: s = 30. Lines are guides to the eye. 

We observe that the spatial decay of the breather is 
exponential for small distances from the center, while it 
becomes algebraic (in fact exactly l/l s ) after a crossover 
at some distance l c (see Fig. 2). 
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FIG. 2. Same as in Fig.l but only for s — 30 in a log-log 
plot. 

Note that l c is s-dependent. Moreover, l c is also de- 
pending on the parameter which selects a given breather 
solution from its one-parameter family (this parameter 
could be the breather frequency, its energy, action or 
something else). In order to understand this result we 
can proceed along the following path. Since the breather 
amplitude decays to zero with increasing distance from 
the center, we can linearize the equations of motion far 
from the breather center, keeping the information that 
we deal with a time-periodic solution with frequency loi, 
fulfilling the nonresonance condition kuj^ ^ £l q (see e.g. 
for details and also for exceptional nonlinear correc- 
tions which are however not important at this stage). 
Assuming Xi(t) = J2k xi(k)e lkuJbt we find linear differ- 
ence equations for xi(k) which do not mix in fc-space. 
The spatial decay of the fc-th amplitude is then given by 
I the lattice Green's function 

[E 

* Ga ^ = / 02 \ d g ' A = k ■ (3) 

The spatial decay of the breather is thus given by the 
convergence properties of the Fourier series, whose coef- 
ficients are given by the rhs in (|^). As is known, nonana- 
lytic functions with discontinuities in the (s— l)-st deriva- 
tive (cf. the integrand on the rhs in (pj)) produce Fourier 
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series which converge algebraically I /I s From that 

follows that at large distances the spatial decay of the 
breather will be algebraic, which is what we found in 
Fig.l. To obtain the exponential decay at small dis- 
tances, let us first slide along the breather family such 
that the breather frequency (or one of its multiples) ap- 
proaches the edge of the phonon band f2 g . Then the 
integrand (|^) will become very large for wave numbers 
close to the band edge which is approached. Applying 
a stationary phase approximation to (|^), i.e. expanding 
the integrand around the band edge we obtain 



G x (l) ' 
for s > 3 and 
Gx(l) 



cos(ql) 



V2-X + C({s - 2)q- 
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cos(ql) 



o v 2 -X + 2Ca(s)q s 



—dq 
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for 1 < s < 3. Standard evaluation of (Q) (closing the 
integration contour in the complex plane by adding a 
half circle with infinite radius and evaluating the residua) 
yields G\( l) ~ e~ %/ ' U2 ~ A/ for s > 3, i.e. exponential de- 
cay [^2|,^3|! On the other side, ((H) yields (closing the 
integration contour in the complex plane by adding a 
quarter circle and returning to zero along the positive 
imaginary axis, and noticing that there are no poles of 
the integrand in the enclosed first quadrant including the 
imaginary axis) G\(l) ~ l/l s for 1< s < 3 @|u|. 

Now we can explain the observed crossover from expo- 
nential to algebraic decay in Fig.l. Indeed, the stationary 
phase approximation for these cases leads to (Q) in the 
limit (i>2 — A) — > 0. This approximation neglects higher 
order terms in the expansion of E s {q) around q = which 
necessarily contain nonanalytic terms. Consequently (^) 
probes (^) for not too large distances (this is counterin- 
tuitive to the assumption that the stationary phase ap- 
proximation is correct for large I [ p"2||T3j ], which it is not). 
Thus we can explain the observed crossover. We can 
even estimate the crossover distance l c using a simple ar- 
gument. A tagged site with index I < l c and I > (the 
center of the breather is located at lb = 0) will experi- 
ence forces from all other sites with index I' according 
to ([!]) . The amplitude of these forces will monotonously 
decay to zero for increasing I' with I' > I, However the 
amplitude of the forces for decreasing /' will be given by 
(I — l')~ s e l '^ l ~ l > for < I' < I (here v is the given expo- 
nent of the spatial decay for \l\ < l c ). Since for negative 
I' the amplitude of these forces will again monotonously 
decay to zero, the worst case is reached when I' = 0. If 
this force acting from the center of the breather on site 
I is comparable to the forces acting on I from its nearest 
neighbours, the exponential decay will be violated. This 
condition yields l~ s e vlc = 1 or 



lnZ e 
I, 



(6) 



This equation has either two solutions or none. For the 
larger (physically relevant) solution we find l c — > oo if 
v I s — > 0, while the smaller one yields 1 in this limit and 
is not of interest. Thus for s > 3 exponential decay is 
reobtained either for large s or for breathers with fre- 
quencies close to the phonon band edge. Since we are 
considering a lattice, the exponential decay part will dis- 
appear if l c *» 1 or smaller. For s — 20 and v — 4.2724 we 
obtain l c = 11.39, and for s — 30 and the same value of v 
the result is l c = 21.56. We miss the observed crossovers 
in Fig.l by just two sites. 

For 1 < s < 3 no exponential decay is observed pro- 
vided the breather frequency is located in the gap be- 
low the phonon band. For breather frequencies above 
the phonon band the dispersion at the upper band edge 
yields always quadratic dependence in q (see above) and 
thus there will be always a crossover from exponential to 
algebraic decay (provided l c > 1). All these results were 
verified by calculating corresponding breather solutions. 

To conclude this part we want to stress that a mod- 
ified interaction 4>2{l) ~ (— l) l /l s will simply exchange 
the notation of upper and lower phonon band edges, and 
the case of acoustic interactions is obtained by letting 

Let us now turn to the question of energy thresholds 
for breathers in the presence of long range interactions. 
There are two lines of argumentation known from the 
short range interaction case |Q . The first one estimates 
the far distance energy of the breather solution in the 
limit when the amplitude of the breather center is small 
and thus its frequency is close to a phonon band edge (the 
only limit where the breather energy can actually become 
small). Using exponential spatial decay the result is that 
the breather energy tends to zero only if ^ or/and 
«4 0, stays finite if — — and v$ ^ or/and 
vq ^ 0, and diverges if V3 = U4 = U5 = vq — and 
Vfi 7^ for some fx > 7 (see Jl^] for details). In the case 
of long range interaction l c tends to infinity in this limit 
for tub > £l q or LUb < £l q and s > 3. Consequently the 
breather energy will have the same qualitative behaviour 
as in the case of short range interaction (the results are 
similar to those obtained in ju) with the tendency that 
the height and the position of the energy minima shift to 
larger values with decreasing values of s). 

However for uib < Q q and 1 < s < 3 no exponen- 
tial decay is observed and the far distance energy of the 
breather is given by ~ A 2 J ^d d r where A is the am- 
plitude of the breather center. This energy will always 
vanish in the limit of zero amplitude. However we are in 
posess of a second line of argumentation for the behaviour 
of the breather energy at small amplitudes. For that we 
consider a finite system of N sites. As was shown in 
fTsfl , the band edge plane waves (BEPW) (which can be 
rigorously defined in the limit of vanishing amplitudes) 
undergo tangent bifurcations, which give birth to dis- 
crete breathers. The amplitude A c of the BEPW at the 



3 



bifurcation point (for nonvanishing cubic and/or quartic 
terms in the Hamiltonian) was calculated in [IS}] : 



O 2 - SI 2 

lL BEPW s '«il 



(7) 



where q\ denotes the wavevector closest to the band edge 
wavevector. Here we consider periodic boundary condi- 
tions and a cutoff in the interaction at one half of the 
system size. This cutoff will induce finite size corrections 
to the dispersion J! 2 for all q except for the band edge 
points. With q% = ^ this correction amounts to 



1 / ?7T V 

A qi =nl(^)-nl(N) = 2C £ _^l-cos(-m)j 



m=JV/2+l 



Evaluation of for s > 1 gives 



AT 
(8) 

(9) 



with 6(s) = 
result for 



"4(1- 

X s v 

and 1 



-cos x)dx. Consequently the correct 



< s < 3 is At 



z{s)/N s 



with 

c ( s ) — Jo J^O- ~ cosa;)dx. The total energy E c ~ YA 2 
in the bifurcation point for 1 < s < 3 is finally given by 



E c ~ N 



(10) 



This has to be contrasted to the case of short range in- 
teractions in one-dimensional systems which can be ob- 
tained from ( |l0| ) by choosing s = 3 and is E c ~ 1/N. We 
thus find, that anomalous dispersion at the band edge 
~ q s ~ l for 1 < s < 3 even further supports the diver- 
gence of the breather energy at small amplitudes, since 
for cubic and quartic anharmoncities in the Hamiltonian, 
for which no divergence in energy is found for short range 
interactions, energy divergence is obtained for long range 
interaction with s < 2. These results confirm studies 
of nonlinear Schrodinger chains with long range interac- 
tions, where s < 2 marks the appearance of two stable 
soliton solutions compared to one for s > 2 |l2| . 

Let us discuss the results. First, we numerically con- 
firm that discrete breathers persist in the case of long 
range interactions, even in the case of anomalous disper- 
sions at the band edge. Secondly, the spatial decay of 
breathers is characterized by a crossover length which 
separates exponential from algebraic decay. Third, we 
show that the existence of energy thresholds for breather 
solutions is supported by long range interactions, and can 
take place when short range interactions (e.g. in one- 
dimensional systems) are not capable of producing these 
thresholds. Thus we can state, that discrete breathers 
appear independent of the lattice dimension and survive 
acoustic and anomalous dispersions. Discrete breathers 
have energy thresholds provided the lattice dimension is 



large enough (typically d > 2) or the interaction is long 
ranged (for d — 1 s < 2). 

Let us speculate on the value of these energy thresholds 
for the lattice dynamics of crystals. As there is no small 
parameter in the system, these threshold energies will be 
comparable to the energy of a vacancy. Consequently 
discrete breathers in two- or three-dimensional crystal 
lattices will be high energy excitations, which could play 
a role close to the melting transition. In contrast for 
one-dimensional systems there are no energy thresholds 
(except for s < 2) and breathers can play a role in the 
dynamics of molecules and similar objects also at low 
temperatures. 



Acknowledgements . 
I wish to thank D. Bonart and J. B. Page for many dis- 
cussions which initiated these studies, C. Baesens and R. 
S. MacKay for sending me their work prior publication 
(and drawing my attention to the helpful reference [ jT5| ) 
and Yu. Gaididei, M. Katsnelson, Yu. Kosevich and O. 
Yevtushenko for helpful discussions. 



[1] S. Aubry. Physica D, 103:201, 1997. 

[2] S. Flach and C. R. Willis. Phys. Rep., 295:181, 1998. 

[3] A. A. Ovchinnikov. Zh. Eksp. Teor. Fiz. / Soviet Physics 

JETP, 57 / 30:263 / 147, 1969 / 1970. 
[4] A. A. Ovchinnikov and H. S. Erikhman. Uspekhi Fiz. 

Nauk (russian), 138:289, 1982. 
[5] S. Aubry, S. Flach, K. Kladko, and E. Olbrich. Phys. 

Rev. Lett, 76:1607, 1996. 
[6] S. Flach and V. Fleurov. J. Phys.: Cond. Mat, 9:7039, 

1997. 

[7] W. Z. Wang, A. R. Bishop, J. T. Gammel, and R. N. 

Silver. Phys. Rev. Lett, 80:3284, 1998. 
[8] R. Livi, M. Spicci, and R. S. MacKay. preprint, 1997. 
[9] S. Flach, K. Kladko, and S. Takeno. Phys. Rev. Lett, 
79:4838, 1997. 
[10] D. Bonart. Physics Letters, A23P201, 1997. 
[11] C. Baesens and R. S. MacKay. Preprint, 1998. 
[12] Yu. B. Gaididei, S. F. Mingaleev, P. L. Christiansen, and 

K. O. Rasmussen. Phys. Rev., E55:6141, 1997. 
[13] S. Flach. Physica D, 113:184, 1998. 

[14] S. Flach, K. Kladko, and R. S. MacKay. Phys. Rev. Lett, 
78:1207, 1997. 

[15] A. Zygmund. Trigonometric Series. Cambridge Univer- 
sity Press, 1968. 

[16] A. A. Maradudin, E. W. Montroll, and G. H. Weiss. The- 
ory of Lattice Dynamics in the Harmonic Approximation. 
Academic Press, New York, 1963. 

[17] J. L. Marin and S. Aubry. Nonlinearity, 9:1501, 1996. 

[18] S. Flach. Physica, D9P223, 1996. 



4 



